Re normalization-group approach to the stochastic Navier— Stokes equation: 

Two-loop approximation 
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The field theoretic renormalization group is applied to the stochastic Navier-Stokes equation that 
describes fully developed fluid turbulence. The complete two-loop calculation of the renormalization 
constant, the j3 function, the fixed point and the ultraviolet correction exponent is performed. 
The Kolmogorov constant and the inertial-range skewness factor, derived to second order of the e 
expansion, are in a good agreement with the experiment. The possibility of the extrapolation of the 
e expansion beyond the threshold where the sweeping effects become important is demonstrated on 
the example of a Galilean-invariant quantity, the equal-time pair correlation function of the velocity 
field. The extension to the d-dimensional case is briefly discussed. 
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I. INTRODUCTION 



One of the oldest open problems in theoretical physics is that of describing fully developed turbulence on the basis 
£^ , of a microscopic model. The latter is usually taken to be the stochastic Navier-Stokes (NS) equation subject to an 
external random force that models the energy injection by the large-scale modes; see, e.g., [fy]|]. The aim of the theory 
is to verify the basic principles of the celebrated Kolmogorov-Obukhov phenomenological theory, study deviations 
from this theory, determine the dependence of various correlation functions on the times, distances, external (integral) 
and internal (viscous) turbulence scales, and derive the corresponding scaling dimensions. Most results of this kind 
were obtained within the framework of numerous semiphcnomcnological models that cannot be considered to be the 
basis for construction of a regular expansion in certain small (at least formal) parameter (l],|^] . 

An important exception is provided by the renormalization group (RG) method that was earlier successfully applied 
in the theory of critical behavior to explain the origin of critical scaling and to calculate universal quantities (critical 
dimensions and scaling functions) in the form of the e expansions jflQ . 

The RG approach to the stochastic NS equation, pioneered in allows one to prove the existence of the infrared 

[~^. ' (IR) scale invariance with exactly known "Kolmogorov" dimensions and the independence of the correlation functions 
of the viscous scale (the second Kolmogorov hypothesis), and calculate a number of representative constants in a 
reasonable agreement with experiment. Detailed review of the RG theory of turbulence and more references can be 
found in |@. 

The standard RG formalism is applicable to the stochastic NS equation if the correlation function of the random 
force is chosen in the form oc fc 4 ~ d_2e , where k is the momentum (wave number), d is the space dimension and the 
£h ' exponent e plays the part analogous to that played by 4 — d in the RG theory of critical behavior. Although the 
results of the RG analysis are reliable and internally consistent for small e, the possibility of their extrapolation to the 
• , real value e = 2 (see below) and thus their relevance for the real fluid turbulence was called in question in a number 
of studies, e.g. 

The crucial role in recent studies of intermittency and anomalous scaling of fluid turbulence was played by a simple 
model of a passive scalar quantity advected by a random Gaussian field, white in time and self-similar in space, 
the so-called Kraichnan's rapid-change model There, for the first time the existence of anomalous scaling was 
established on the basis of a microscopic model ]l7| and the corresponding anomalous exponents were calculated 
within controlled approximations Jl^Jl9[| and numerical experiments [ po[ . Detailed review of the recent theoretical 
research on the passive scalar problem and more references can be found in pl[ . 

The field theoretic RG and the operator product expansion (OPE) were applied to the rapid-change model in [p2[ . 
The RG allows one to construct a systematic perturbation expansion for the anomalous exponents, analogous to the 
e expansion in models of critical behavior, and to calculate the exponents to the second p3] and third |23| orders. 
For passively advected vector fields the exponents for higher-order correlation functions were derived by means of the 
RG techniques to the leading order in e in Refs. pifl . It was shown that the knowledge of three terms allows one to 
obtain reasonable predictions for finite e ~ 1; even the plain e expansion captures some subtle qualitative features of 
the anomalous exponents established in numerical experiments p3| |. 

In contrast to standard </> 4 model of critical behavior J3|,E| , where the critical exponents are known up to the order 
e 5 (five-loop approximation), and to the rapid-change model, where the anomalous exponents are known up to e 3 , 
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all the calculations in the RG approach to the stochastic NS equation have been confined with the simplest one-loop 
approximation. The reason for this distinction is twofold. First, any multiloop calculation for this dynamical model 
is a demanding job: one can say that the two-loop calculation for the stochastic NS equation is as cumbersome 
as the four-loop calculation for the conventional 4 model. Second, the critical dimensions for the most important 
physical quantities (velocity and its powers, frequency, energy dissipation rate and so on) are given by the one-loop 
approximation exactly (the corresponding e series terminate at first-order terms) and the higher-order calculations 
for them are not needed. 

However, the e series for other important quantities do not terminate and the calculation of the higher-order terms 
for them is of great interest. In this paper, we present the results of the two-loop calculation for a number of such 
quantities: the j3 function, the coordinate of the RG fixed point, the ultraviolet (UV) correction exponent u>, the 
Kolmogorov constant Ck and the inertial-range skewness factor S. The knowledge of higher-order terms is also 
important to judge about the validity and convergence properties of the e expansions in the models of turbulence on 
the whole. 

The plan of the paper is as follows. In Sec. [n| we introduce the stochastic NS equation paying special attention to the 
choice of the random force correlator and the physical value of the parameter e. In Sec. Ill we recall the field theoretic 
formulation of the model, diagrammatic technique, renormalization and RG equations. Since the "ideology" of the 
RG and OPE approach to the stochastic NS model is explained in Refs. ||,[l(],^5|] in detail, here we confine ourselves 
to only the necessary information. In Sec. IV we present in detail the two-loop calculation of the renormalization 
constant, RG functions, coordinate of the fixed point and the UV correction exponent. We explain some tricks (e.g. 
an appropriate choice of the IR cutoff) that essentially simplify the calculation and make it feasible. 

In Sec. we calculate two orders of the e expansion for the pair correlation function of the velocity field (this 



accuracy is consistent with the two- loop calculation of the RG functions in Sec. IV). They are used later in Sec. VII 
in the calculation of the skewness factor and the Kolmogorov constant. 

In Sec. VI we discuss on the example of the pair correlation function the possibility of the extrapolation of the e 
expansion to finite values of e. We explain an important distinction that exists between the RG technique and the 
traditional approach based on diagrammatic self-consistency equations and argue that, contrary to what is sometimes 
claimed, the RG is applicable beyond the threshold where the so-called sweeping effects become important, but one 
should co mbin e the RG with the OPE and go beyond the plain e expansions. 

In Sec. VII we perform the two-loop calculation of the Kolmogorov constant Ck and the inertial-range skewness 
factor S. The new point is not only the inclusion of the second-order correction, but also the derivation of Ck 
through an universal (in the sense of the theory of critical behavior) quantity. This allows us to circumvent the main 
drawback of earlier calculations of Ck- the intrinsic ambiguities in the corresponding e expansions. To the best of 
our knowledge, the inertial-range skewness factor S has never been calculated earlier within the RG approach. 

The experience on the RG theory of critical behavior shows that for dynamical models the higher-order corrections 
are not small, and in order to obtain reasonable predictions, say, for critical exponents, one should augment the 
plain e expansions by additional information obtained from the instanton calculus and Pade-Borel summations; the 
situation for the amplitudes in scaling laws is even worse (see e.g. ||). Surprisingly enough, in our case already 
the first-order results are in a reasonable agreement with experiment, which is not destroyed by the inclusion of the 
second-order corrections. It turns out that the experimental value of Ck and S lie in between of the two consecutive 
approximations, like for the exactly solvable Heisenberg model p6| . 

The main conclusion of the paper is that the renormalization group and the e expansion allow one to derive a 
number of important characteristics of the real fluid turbulence beyond the simplest first-order approximations and 
in a good agreement with the experiment. 



II. STOCHASTIC NAVIER— STOKES EQUATION AND THE CHOICE OF THE RANDOM FORCE 

As the microscopic model of the fully developed, homogeneous, isotropic turbulence of an incompressible viscous 
fluid one usually takes the stochastic NS equation with a random driving force 

V t <Pi = vo&<Pi - diP + ft, V t = d t + (ipd). (2.1) 

Here ipi is the transverse (due to the incompressibility) vector velocity field, V and /, are the pressure and the 
transverse random force per unit mass (all these quantities depend on x = {i,x}), isq is the kinematical viscosity 



coefficient, d 2 is the Laplace operator and Vt is the Lagrangian derivative. The problem (2.1) is studied on the entire 
t axis and is augmented by the retardation condition and the condition that ipi vanishes for t — > — oo. We assume for 
/ a Gaussian distribution with zero average and correlator 
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{h{x)f 3 {x')) = S -^-p- J dkP tJ (k)d f (k) exp [ik(x-x')], (2.2) 

where Py(k) = <5y — kikj/k 2 is the transverse projector, df(k) is some function of k = |k| and model parameters, and 
d is the dimension of the x space. The time decorrelation of the random force guarantees Galilean invariance of the 



full stochastic problem (2.1), fl2.2|). 

Let us specify the form of the function df{k) in the correlator ( [2.2| ) used in the RG theory of turbulence; more detailed 
discussion is given in Refs. Physically, the random force models the injection of energy to the system owing 

to interaction with the large-scale eddies. Idealized injection by infinitely large eddies corresponds to df(k) oc <5(k), 
more precisely, 

d f (k) = 2(2ir) d £S(k)/(d- 1), (2.3) 

where £ is the average power of the injection (equal to the average dissipation rate) and the amplitude factor comes 
from the exact relation 

^ = ^4r / dkd f (h). (2.4) 



2(27r) d 

On the other hand, for the use of the standard RG technique it is important that the function df(k) have a power-law 
behavior at large k. This condition is satisfied if df(k) is chosen in the form || 

d f (k)=D k±- d - 2e , (2.5) 

where Dq > is the amplitude factor and e > is the exponent with the physical value e = 2 (see below). 
Let us recall the well-known power-law representation of the <i-dimensional S function, 

6(k) = lim j^Jdx (Ax) 2 - 4 exp[i(k • x)] = Sj 1 ^ Urn [(4 - 2e)(fc/A) 4 " 2e ] , (2.6) 

with some ultraviolet (UV) momentum scale A. Here and below we denote 

S d = 2ir d / 2 /T(d/2), S d = S d /(2n) d , (2.7) 
where S d is the surface area of the u nit s phere in c?-dimensional space and T(- • ■) is Euler's Gamma function. It then 



follows that for e — > 2, the function (|2.5|) turns to the ideal injection (2.3) if the amplitude Dq is related to £ as 



^ 4(2_ £)A 2£ 4 - for £ — > 2. (2.8) 



A more realistic model, used e.g. in ||[25|], is 

d f (k) = D k 4 - d ~ 2e h(m/k), h(0) = 1, (2.9) 

where m — 1/L is the reciprocal of the integral turbulence scale L and h(m/k) is some well-behaved function that 
provides the IR regularization. Its specific form can be chosen to simplify the practical calculation in higher orders; 



see Sec. [V 



In the RG approach to the problem (2.1), ( |2. 2]^ ( |2.9| ) the exponent e plays the part analogous to that played by 
4 — d in Wilson's theory of critical phenomena [[Hj]. All the results of the RG analysis are reliable and internally 
consistent for small e: the possibility of their extrapolation to the real value e = 2 was called in question in a number 
of studies, e.g. ||ll| - |l5|. We shall return to this im por tant issue later, and here we only note that the replacement of 
the ideal injection ( |2.3| ) by the power-law model (Ej) followed by the use of the e expansion was confirmed by the 
example of the exactly solvable Heisenberg model pa] . 

III. FIELD THEORETIC FORMULATION AND RENORMALIZATION OF THE MODEL 

Detailed exposition of the RG theory of turbulence and the bibliography can be found in ||,[l0| ; below we restrict 
ourselves to only the necessary information. 
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Stochastic problem (2.1), is equivalent to the field theoretic model of the doubled set of transverse vector 

fields $ = {ip, ip'} with action functional 



= if' D ftp' /2 + <f/[-8t<p + v d 2 v - (<pd)ip], 



(3.1) 



where Df is the random force correlator (2.2) and the required integrations over x = {£,x} and summations over the 
vector indices are understood. 

The model ( |3.l| ) corresponds to a standard diagrammatic perturbation theory with bare propagators 

( w ')o = (pV)o = + wbfc 3 )- 1 , (^V')o = 0, ( W )o = rf/(fc)/(^ 2 + "ok 4 ) (3.2) 
in the frequency-momentum (u>— k) representation or 



(^(ty(* / ))o = fl(*-0® c P 



-i>ok 2 (t 



t 1 )} , (tp'{t)tp'(f)) = 0, ( W )o = exp [-^ fc 2 |f - t'|] (3.3) 



in the time-momentum (t-k) representation; the common factor Py(k) in (|3.2|), (3.3) is understood. The interaction 



in (|3.l| ) corresponds to the triple vertex —ip'(ipd)ip — ip' i Vij S (fjip s /2 with vertex factor Vy s = i(kjSi S + fc s %), where k 
is the momentum of the field 99'. The part of the cou pling constant (expansion parameter in the ordinary perturbation 
theory) is play ed b y go = Dq/i/q with Dq from (2.9). 

The model (3.1) is logarithmic (the coupling constant go is dimensionless) at e = 0, and the UV divergences have 
the form of the poles in e in the correlation functions of the fields $ = {<p,f'}. Superficial UV divergences, whose 
removal requires counterterms, are present only in the 1-irreduciblc function (tp'tp), and the corresponding counterterm 
reduces to the form ip'd 2 ip. In the special case d = 2 a new UV divergence appears in the 1-irreducible function (tp'(p'}. 
We shall study this important case in a separate paper, and for now we always assume d > 2. Then for the complete 
elimination of the UV divergences it is sufficient to perform the multiplicative renormalization of the parameters i/q 
and go with the only independent renormalization constant Z v : 



vo 



vZ,, 



9o 



2e v 



(D 



go^o 



(3.4) 



Here /1 is the reference mass in the minimal subtraction (MS) scheme, which we always use in what follows, g and v 
are renormalized analogs of the bare parameters go and vq, and Z = Z(g, s, d) are the renormalization constants. No 
renormalization of the fields and the "mass" mo = m is needed, i.e., Z$ = 1 for all <£> and Z m = 1. 
In the MS scheme the renormalization constants have the form "1 + only poles in e," in particular, 



Z u = l + ^a k {g)e fe = l + ^ 5 n ^a„ fe e k , 



fe=l 71=1 fe=l 

where the coefficients a n k depend only on d. The one-loop result 



(3.5) 



an 



-(d-l)S d /8(d + 2), 



(3.6) 



with Sd from (2.7), was presented in ||. 

Since the fields are not renormalized, their renormalized correlation functions W R coincide with their unrenormalized 
analogs W = ($...$); the only difference is in the choice of variables and in the form of perturbation theory (in 
g instead of <7o) : W (g, v, fi, m, . . .) = W(go, vo, m o, ■ ■ ■)■ Here the dots stand for other arguments like coordinates, 
times, momenta and so on. We use 2? M to denote the differential operator fid^ for fixed bare parameters go, vo,mo 
and operate on both sides of that relation with it. This gives the basic differential RG equation: 

V RG W R (g,v,n,m,...) = 0, V RG = X> M + /3(g)d g - >y v (g)T> v , (3.7) 



where T>rg is the operation T>^ expressed in renormalized variables, T> x = xd x for any variable x, and the RG functions 
(the anomalous dimension 7^ and the f3 function) are defined as 



lv {g) = V^\nZ v = -2V g a x {g), 0(g,e) = V^g = g [-2e + 3j v (g)} 



(3.8) 



with ai(g) from (3.5). In the MS scheme only the residues at the first-order poles in s, that is, only the coefficients 
afci, contribute to the RG functions owing to the UV finiteness of the latter; this explains the last relation for 7„. 
The relation between 8 an d 7^ results fro m the definitions and the last relation in (3.4). 

From the relations (3.8) using Eq. (3.6) we find the first-order expressions for the RG functions: 
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lv {g) = (d- l)S d g/A{d + 2) + 0(g 2 ), [3(g, e) = g [-2s + 3(d - l)S d g/4(d + 2)] + 0(g 3 ). 



(3.9) 



From (3.9) it follows that an IR-attractive fixed point 

g* = 8(d + 2)e / 3(d - + 0(e 2 ), /%*) 



0, uj = f3'{g*) = 2e + O{e 2 )>0 



(3.10) 



of the RG equation ( [Oh ex ists in the physical region g > for e > 0. The value of J v (g) at the fixed point is found 
exactly using relations (p.8[): 



(3.11) 



without corrections of order e 2 , e 3 , and so on. 

For definiteness, consider the solution of the RG equation (3.7) on the example of a correlation function W R 
that involves n fields tp and depends on a single coordinate difference r = |r|; the extension to the general case 
is straightforward. We shall omit the superscript R in what follows. In renormalized variables, dimensionality 
considerations give: 



W = {v/r) n R(s,u,g), 



txr, 



where R(- • •) is a function of completely dimensionless arguments (the dependence on d and e is understood). 
From the RG equation the identical representation follows, 

W={u/r) n R{s,u,g) = (9/r) n R(l,u,g), 



(3.12) 



(3.13) 



where the invariant variables e = e(s, e) satisfy the differential equation T>jiqE = with the operator T>rg from (3.7) 
and the normalization conditions e = e at s — 1 (we have used e = {v, g,m} to denote the full set of renormalized 
parameters). The identity u = u is a consequence of the absence of T> m in the operator T>rg owing to the fact that m 
is not renormalized. Equation ( 3.13| ) is valid because both sides of it satisfy the RG equation and coincide for s = 1 
owing to the normalization of the invariant variables. 

The relation between the bare and invariant variables has the form 



v = vZ v (g), 



9o = gr 



- 2e Z g (g) 



(3.14) 



with from (3.4). Equation (3.14) determines implicitly the invariant variables as functions of the bare parameters; 
it is valid because both sides of it satisfy the RG equation, and because Eq. (3.14) at s = [ir = 1 coincides with (3.4) 
owing to the normalization conditions. 

It is well known th at th e behavior of the invariant charge g at large \xr is governed by the IR stable fixed point : 
g — > 5* with from (|3.10|). Then the large-/ir behavior of the invariant viscosity v is found explicitly from Eq. (3.14) 
and the last relation in (|3.4| ): 



Then for s — > oo and any fixed u 



v = [D Q r 2 yg] 
mr we obtain 



1/3 



[D r 2e /g< 



,1/3 



W=(Do/g*) n/3 r- nA «f(u), 



A.. 



l-2e/3, f(u) = R(l,u,g,). 



(3.15) 



From Eq. (3.15) it follows that in the IR range (large \ir and arbitrary mr) the parameters go and vq enter into the 
correlation functions only in the form of the combination Dq = g^ v^, a fact first established in Ref. 0. For e = 2 this 
proves the Second Kolmogorov hypothesis: owing to the relation ( |2.8| ), the correlation functions depend on £ but do 
not depend on the visc osity coefficient vq and the UV scale A ~ So; see also the discussion in |Tc|^^| . 

Representation ( 3.15| ) for any scaling function f(u) describes the behavior of the correlation functions for s = fir 3> 1 
and any fixed value of u = mr; the inertial range correspo nds to the additional condition u <C 1 . T he form of 
the function f(mr) is not determined by the RG equation ([Tt]). Calculating the function R in (3.12) within the 
renormalized perturbation theory, R = ^2^— g n Rn, substituting g — > g* and expanding g* and R n in e, one obtains 
the e expansion for the scaling function: 



/(«) 



n=0 



S n fn(u). 



(3.16) 



Although the coefficients /„ in our model are finite at u — 0, this does not prove the finiteness of f(u) beyond the 
e expansion: one can show that for any arbitrarily small value of e there are diagrams that diverge at m oc u — > 0. 
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As a result, the coefficients f n contain IR singularities of the form u p \n q u, these "large IR logarithms" compensate 
for the sm allnes s of s, and the actual expansion parameter appears to be elnu rather than e itself. Thus the plain 
expansion ( 3.1 6| ) is not suitable for the analysis of the small-it beh avior of /(it). 



The formal statement of the problem is to sum up the expansion (3T£) assuming that e is small with the additional 
condition that elnu ~ 1. By analogy with the theory of critical behavior the desired solution can be obtained 

using the well-known operator product expansion (OPE); see |9|, |l0| , p5| . The OPE shows that for small u the function 
f(u) has the form 

/(u) = ]Tc F (u) U A -, (3.17) 

F 

where Cp are coefficients analytical in u 2 , the summation runs over all local composite operators F(x) consistent with 
the symmetries of the model and the left-hand side, and are their scaling dimensions calculated as series in e. 



If all Ap are non-negative (as in models of critical behavior), expression (3.17) is finite at u — and the e expansion 
can be used for the calculation of /(0). The feature specific to statistical models of fluid turbulence is the existence 
of composite operators with negative dimensions Ap. Such operators, termed "dangerous" in [p pio|j25[ |, give rise to 
strong IR singularities in the correlation functions and to their divergence at m cx u — > 0. 

The problem is that dangerous operators in our model are absent in the e expansions and can appear only at finite 
values of e. This means that they can be reliably identified only if their dimensions are derived exactly with the aid of 
Schwinger equations or Ward identities that express Galilean symmetry. Due to the nonlinear nature of the problem, 
they enter the corresponding OPE's as infinite families whose spectra of dimensions are not bounded from belo w, and 



in order to find the small-mr behavior one has to sum up all their contributions in the representation ( 3.17 ). The 
needed summation of the most singular contributions, related to the operators of the form (p n with known dimensions, 
was performed in |25|] using the so-called infrared perturbation theory for the case of the different-time pair correlation 
functions; see also P,^0|- It has revealed their strong dependence on m, which physically can be explained by the 
infamous "sweeping effects." This demonstrates that, contrary to the existing opinion |ll|-|l5|, the sweeping effects 
can be properly described within the RG approach, but one should combine the RG with the OPE and go beyond 
the plain e expansions. 

In Sees. M-WTl we shall be interested in the Galilean invariant objects like the equal-time structure functions. For 



such objects, the singular contributions mentioned above drop out from the representation (3.17) in agreement with 



the fact that Galilean invariant quantities are not affected by the sweeping. Only the dimensions of Galilean invariant 



operators can appear in (3.17), and the singular behavior can be related to invariant operators with Ap < 0. No such 
operator, however, has been presented in any study we know of; see Refs. ||^7],^8| . Although the problem shoul d be 



considered open, this fact suggests that Galilean invariant objects may have a finite limit at u — > 0. Below in Sec. VII 
we shall calculate the Kolmogorov constant within the assumption that the second-order structure function has a 
finite limit at u — ► for the physical value e = 2. Such an assumption can be justified by the real experiments, which 
indicate that the corresponding exponent is hardly distinguishable from the Kolmogorov value Q , or by the example 
of the Kraichnan model, in which the second-order function is not anomalous fll6[. 



IV. TWO-LOOP CALCULATION OF THE RENORMALIZATION CONSTANT, RG FUNCTIONS, 
FIXED POINT AND THE UV CORRECTION EXPONENT 



Let us turn to the calculation of the renormalization constant Z v in Eq. (3.5) with the accuracy 0(g 2 ) (two-loop 
approximation) . 

The renormalization constants in the M S sc heme are independent of the value of the "mass" m and the specific form 
of the f unct ion h(m/k) in the correlator (2.9), and the practical calculations are usually performed in the "massless" 
model (2.5) with m = an d h (m/k) = 1. It turns out, however, that for the multiloop calculations it is more 
convenient to use the model ( [2.9| ) with m/0 and h(m/k) = 6{k — m). We shall determine the constant Z v from the 
requirement that the 1-irreducible correlation function (ip'<p) at zero frequency (w = 0) and i n th e limit p — > be UV 
finite (that is, be finite at e — > 0) when expressed in renormalized variables using relations ([O]). More precisely, we 
shall consider the scalar dimensionless quantity 



lim 



(^i)i_i r (p; w = °) 
vp 2 {d — 1) 



which depends only on g and mj \i and in the loopless (tree) app roximation equals r^ -* 
is provided by the sharp cutoff h(m/k) = 8(k — m) in Eq. (|2.9D. 



(4.1) 



-Z v . The IR regularization 



G 



Let us illustrate this scheme by revisiting the one-loop calculation. To first order in g, the quantity (4.1) is 
determined by the only diagram 




(4.2) 



Here the propagators (3.2) are represented by the lines, the slashed ends correspond to the field if', the ends without 
a slash correspond to <p. The triple vertices correspond to the vertex factor Vij s (see the text below Eq. (3.3)), and 
the contraction with the transverse projector is implied. Integrating over the frequency and taking the limit p — > 
gives 



r 



(i) 



2e 7-1 



rfk 



4(d- l)(27r) d J fc d+2e 



h(m/k) [d - 3 + (9 - d)z 2 - 6z 4 



(4.3) 



where z = (n • k)/fc is the cosine between k and n = p/ p, th e direction of the external momentum. The quantity T 
does not depend on the direction n and we can average (4.3) with respect to it using the relations 



(z 2n ) = 



(2n- 1)!! 



y 2n+l\ _ 



= 0; 



d(d+2)...(d + 2n- 2) v 
the brackets denote the averaging over the unit sphere in d dimensions. This gives 



m _ g^S d {d-l) 
4(d + 2)Z 2 



dk 



(d-l)^/m)^gS d 
8(d+2)sZZ 



with S d from (2.7). In the expansion in the number of loops, 

r = -z„ + r (1) + r (2) 



(4.4) 



(4.5) 



(4.6) 



the n-loop term is the quantity of order 0(g n ), so that in the O(g) approximation in (4.E) it is sufficient to keep 
the ter m I^ 1 ) (and put Z v = 1 in it) and the linear in g contribution in the first term. The residue at the pole in e 
in ( f4.5| ) is found by the simple substitution (/i/m) 2e — > 1. From the requirement that the function T(g,m/ fi) be UV 
finite (that is, finit e at e — > 0) we thus find an = — (d — l)Sd/8(d + 2) for the coefficient an in agreement with the 
well-known result (3.6). 

The two-loop contribution in (4.6) is given by the sum of the following eight diagrams (normalization according 
to (O) is implied): 





r 



(2) 




(4.7) 



After two simple integrations over the frequencies, these diagrams are represented as integrals over two independent 
momenta ki, k2 that can be assigned to the two ((f(p)o lines, while the external moment um p flows only via the 
(ip(p')o lines (such choice of the integration contours is possible for any of the diagrams in (4.7)). The way in which 
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the momenta ki and k 2 are assigned to the (f<p)o lines is explicitly shown in (4.7); additional symmetrization in ki, 

(2) . . (2) 

k2 is understood for the diagrams r> with i = 5, . . . , 8. The diagram T\ is automatically symmetric. 

After the contractions of the tensor indices and the limit transition p — » have been performed, the resulting 
integrands take on the form of polynomials in the cosines z\ = (n- k/) /fcx and z 2 = (n • k 2 )/fc 2 . The integrands are 
further simplified after averaging over the direction n = p/p using (4.4), which is possible owing to the independence 



(2 

of the quantities 1^ ; of n. In the variables ki, k 2 and z = (ki • k 2 )/fcifc 2 they can be written in the form 



r (Z> — n l < A iA £ 7 
1 i — 9 ^df 1 A > 



dki 



l+2e 



dk 2 



k. 



l+2e 



dzf i (z,k 1 /k 2 ). 



(4.8) 



The lower limits in the integrals over k\, fc 2 are due to the choice h(m/k) = 9(k — m). For convenience reasons, here 
and in analogous formulas below we restrict the integration over z to the half- interval [0, 1] instead of the full interval 
[— 1, 1] and simultaneously replace the integrands by their doubled even-in- z parts. The e xpr essions for the functions 
fi are given in the Appendix. Within our accuracy, it is necessary to replace Z v — > 1 in (^^) . 
In terms of the dimensionless variables Ki = ki/m one can write 



if' = g*S-jMm)«M£ 



with 



Ms) 



l+2e 



dn 2 

,A+2e 



dzfi{z, ki/kz). 



We are interested in the coefficients a* , bi in the pole part of the quantities Ai (e) : 

£ bi 



Me) 



0(e). 



(4.9) 



(4.10) 



(4.11) 



(2) 

The subdiagram in T\ ' does not contain an UV divergence, so that fi(z, 0) = fi(z, 00) = 0, the integrals over ki 

and fc 2 are separately finite, and the divergence at e — > comes only from the region where k\, k 2 tend to infinity 

(2) (2) 

simultaneously. As a result, the second-order pole is absent: 04 = 0. The subdiagrams in T y 5 . . .,T y 8 ' contain UV 

/ 2 \ 

divergences, but in the sum Iq = 



i=5 



(2) 

T\ ' they cancel each other (a consequence of the absence of divergence in 
the 1-irreducible function (tp'tptp)), so that a = J2i=5 a i — 0- 

For the diagrams Ft with i = 1, 2, 3 o ne still has fi(z, 0) = 0, but the second-order poles exist due to the relations 
fi(z, 00) = const ^= 0. Expressions (4.10) can be simplified using the identity 



V m Tf ] =-Aeg 2 S 2 d (ix/m)^ Me) 



(4.12) 



that follows from (4.9). The operation T> m = md/dm red uces the nu mbe r of integrations in (4.1C) and explicitly 
isolates the pole factor e -1 : calculating the left-hand side of ( 4.12 ) using (4.£) and changing to dimensionless variables 
again, one obtains 



Me) 



1 

47 



dn 

IT+2i 



dz [fi(z, k) + fi(z, 1/k)] . 



For % = 0, 4 the integral in (4.13) is finite at e — and determines the residue at the first-order pole: 



0, 



°° dK 
K 



dz[fi(z,K) + Mz,1/k)], i = 0,4. 



(4.13) 



(4.14) 



For the quantities Aj( e) wi th i = 1, 2, 3 the residue at the second-order pole is obtained if the value of the function 
[fi(z, k) + fi(z, 1/k)] in (4.13) is substituted by its limit value fi(z, oo) + /i(z, 0) = fi(z, 00) [we recall that fi(z, 0) = 0]; 
then the integration over k becomes simple and gives 



1 



The remaining integral with the replacement fi(z,K 
residue at the first-order pole: 



dzfi(z, 00), i = 1,2, 3. 

[fi(z, k) — /i(z, 00)] is finite at e 



(4.15) 

and determines the 
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h = \ / — / dz [fi(z, n) - fi 

4 A K Jo 



(z, oo) +fi(z, 1/k)], i = l,2,3. 



(4.16) 



According to (|4.9|), (4.11), the summed values 



(4.17) 



i=0 



i=0 



of the coefficients a^, &j f rom ( 4.14 )— (4.16) determine the needed residues for r' 2 ) = J2i=o M • Analysis shows that 
all integrations over z in (p^H4re|) can be performed in terms of elementary functions for any d, which gives simple 
expressions for a,. For their summed value a in (4.17) this gives: 



(d~l) 2 
64(d + 2) 2 



<r 2 « 2 
15 d a u 



(4.18) 



with from Eq. (2.7) and an from (3.6). 

It turns out, however, that the most convenient way to find the residues bi is the direct numerical calculation of 
the double integrals in (4.14), (4.16), which can be performed for any given value of d. Below we give the results 
of such calculation for the most interesting case d = 3; it is convenient to represent them in terms of the quantities 
b[ ee 10 3 • b t and b 1 ee 10 3 • b = Eto b 'i- 



b' 



-1.34, b\ = 0.16, b' 2 = 1.21, b' 3 = 0.006, b' A = -4.17, b' 



-4.13. 



(4.19) 



Now let us writ e down the condition of the UV finiteness of the right-hand s ide of Eq. (4.6) in the second order 
in g. In Eq. (4.8) it was sufficient to replace Z v — > 1, while in the expression (4.5) for rW one should keep in the 
expansion Z~ = 1 — 2ang/s + . . . the term linear in g. Then using Eq. (3.5) the condition of the UV finiteness is 
written in the form 



a 22 e 



a 2 ie x ) - 2a\ x £ 2 



1 



2eln-^ 



S d ae 



1 + 4e In — 

m 



Sjbs' 1 = 0. 



(4.20) 



From Eq. ( 4. 18[ ) it is easily seen that the terms with ln(/i/m) in ( 4.2Cf) cancel each other (a c onsequence of the 
renormalizability of the model in the MS scheme). For the sought coefficients a 2 2 and a 2 i in 

(|J) fr om E q s - ( EH > 

and ( 4.18 ) we thus obtain 



a-22 



021 



sib 



(4.21) 



with an from Eq. ( p.6[) and b from ( 4.17 ). This completes the calculation of the renormalization constant Z v in the 
two-loop approximation. 

The knowledg e of the renormalization constant Z v to order 0(g 2 ) allows for the calculation of the anomalous 
dimension 7„ in (3.8) to order 0(g 2 ): 



7„(ff) = -2gd g ai(g) = -2 (a n g + 2a 21 g 2 ) + 0{g 3 ), 



(4.22) 



then the fi function is found from ( p.8[ ) with the proper accuracy. The coordinate of the fixed point in ( 3.10) is then 
determined to order 0(e 2 ). From Eqs. (4.21) and ( 4.22| ) we obtain: 

g„S d = ae{\ + \e) + 0(e 3 ), a = -S d /3a u = 8(d + 2)/3(d - 1), A ee 2a 2 i/3a 2 1 = 2b[8(d+ 2)/{d — l)] 2 /3 (4.23) 

with the coefficient b from (4.17) and Sd from fl2.7| ). In particular, for d = 3 using Eq. ( 4 . 1 9| ) we obtain: 

g ^ = (407r 2 e/3)(l + Ae) + O(£ 3 ), A ~ -1.101. (4.24) 



The case of general d dimensions is also of interest and will be discussed in Sec. VII. Results for the coefficient A for 
different values of d are given in Table Q a long with other quantities that will appear later on. 

Th e corr ection expo nent ui = f3'(g*) in ( J3. 10 ) , determined by the slope of the j3 function at the fixed point, is found 
from ( 4.22] ) and ( [3.10] ) to second order of the e expansion: 



2e(l-Ae)+0(e 3 ). 



(4.25) 



From Table | it follows that, for any d, the two-loop correction makes the value of u> "more positive," that is, it 
enhances the IR stability of the fixed point. 
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V. PAIR CORRELATION FUNCTION: e EXPANSION 



In this section we shall discuss the equal-time pair correlation function of the velocity field in the momentum 
representation, 



G ij (p) = P ij (p)G(p). 



The analogs of representations ( 3.12 ) and ( 3.15 ) for G(p) have the forms 



G(p) 



= gv 2 p~ d + 2 R(s,u,g) ~ I)f 9 yV^i?(l, M ,), 



= p/M, u = m/p, 



(5.1) 



(5.2) 



where the second relation holds in the IR asymptotic region and one factor g is explicitly isolated such that the 
expansion in g of the dimensionless function R(l,u,g) begins with 0(g ). 

Our aim is twofold. Below in this section we shall calculate two orders (1 and e) of the e expansion for the scaling 
function R(l, u, g*) at u = (this accuracy is consistent with the two-loop calculation of the RG functions in Sec. IV). 
The coefficients of the e expansions in our model have a finite limit at u = (see the end of Sec. Ill) so we can 
perform this calculation in the "massless" model (2J5) with m — and h(m /k) — 1, which is always assumed below 
unless stated to be otherwise. These results will be used later in Sec. |VII| in the calculation of the skewness factor 



and the Kolmogorov constant. In Sec. VI, we shall use the example of the pair correlator to discuss the existence of 
the limit u —* for finite £, the possibility of the extrapolation of the e expansion beyond the threshold where the 
sweeping effects become important (see the remark in the end of Sec. Ill), and the relation between the RG approach 
and other approaches to the stochastic NS equation. 

In order to calculate two terms of the s expansion of the amplitude i?(l,0,<7*) we need to find the coefficients c\ 
and C2 in the expansion 



R(l,0,g) = G(p)/ {gv 2 P - d+2 ) |p =/1) „=o = ci + gS d c 2 + Q(g 2 ) 



(5.3) 



(the factor Sd from Eq. ( pTjj ) is isolated in the second term for convenience). We shall see below that c x — 1/2, while 
the coefficient C2 = 02(e) should be found to the leading order in e, that is, at e = 0. 

With the needed accuracy of 0(g 2 ), the function G(p) is given by the sum of the loopless (tree) and one-loop 
diagrams. The first (tree) contribution to the correlator (5.1) is simply given by the bare correlator (ff)o from 
( |3.3[ ) at t = t' and expr esse d in renormalized variables using Eq. (3.4). The corresponding contribution Rq to the 
dimensionless quantity ( |5.3| ) has the form 

R = Z~ l /2 = [1 - ang/e + 0(g 2 )} /2 (5.4) 

with the coefficient an from ( |3.6|) . This gives the exact result c\ = 1/2 for the first coefficient in (|5 
The one-loop contribution to the function G(p) is given by the sum of the following diagrams 




(5.5) 



the diagr amm atic notation was explained below Eq. (4J2). We shall denote their contributions to the dimensionless 
function (5.3) as ill, 2, 3 (from the left to the right). They have the forms 



f dn 



.2— d— 2s 2— d— 2e 



^2-d-2e 



J^)dU(K,q)Q 2 (K)K 2 

1 U(K,q)Q 2 (c L ) q 2 - d - 2£ , 



8(d-l) J (2tt) 
g f dK 



Ri 
R2 

Rs ~ 8(d-l) J (2ttY 

where we included the symmetry coefficient 1/2 for R\ and denoted 

Qi(k) = d - 1 - 2d(n ■ k) + 2(d - 2)k 2 + 4(n ■ k) 2 , 
Q 2 (k) = 1 - d+ 2(d - l)(n ■ k) - (d - i)n 2 - 2 K 2 (n • k), 
U(k, q) = [k 2 - (n • K) 2 ]/K 2 q 2 (l + k 2 + q 2 ). 



(5.6a) 
(5.6b) 
(5.6c) 



(5.7) 
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In Eqs. (5.6) and (5.7) we use only dimensionless variables, namely, the momenta divided by the modulus p of the 
external momentum: n = p/p is the direction of the external momentum, k = k/p is the momentum flowing via the 
{(p<p')o line for i?3, the (<pip)o line for R2 and any one of the two equivalent (<p<p)o lines for R±, and q = n k is the 
mom entum flowing via the remaining line. In order to isolate the scalar factor G(p) we contracted the full expression 
(5.1) with the transverse projector Py(p) and divided the result by its trace; this explains the origin of the factor 
Pu(p) = (d — 1) in the denominator. 

It is easily checked that Qi(k) = Qi(q) and n 2 ~ (n • k) 2 = q 2 — (n • q) 2 . From the last relation it follows that 
U(n,q) = U(q, n) and therefore R2 — -R3. Another consequence of t he last relation is that U(k, q) remains finite 
at k — ► and q — > 0, so that the IR convergence of the integrals ( |5.6| ) at k — > and q — > are determined by the 
factors n 2 ~ d ~ 2£ and q 2 ~ d ~ 2e , respectively. The IR divergences in all three integrals (5.6) aris e for e > 1 , but in their 
sum they cancel each other: the singularities at k — (when q 2 = 1) in the integrands (5.6a) and (5.6b) cancel each 
other due to the relation 



<2i(kOU=o = (d - 1) = -Q2(«0U=o, 
while the singularities at q = (when k 2 = 1) in the integrands ( |5.6aD and (5.6c) cancel out due to the relation 

Qi(/c)| g =o = Qi(q)| q =o = (d-l) = -<2 2 (q)|„=o, 

where we have used the identity Q\(k) = Qi(q). Thus the IR divergences that appear in the individual diagrams 
-Ri,2,3 at e — 1 cancel each other in their sum, and the latter becomes IR divergent only for e = 2. These general 
facts will be illus trated by explicit formulas for the limit d — ► 00 in Sec. fvj| . 



Formulas (5.6) are written for the "massless" model (2.5) that is, h = 1 in Eq. (2.9). In the "ma ssi ve" m odel 
with h =^ 1 additional factors h(u / n)h(u / q) , h(u / K)h[u) , and h(u/q)h(u) appear in the integrands ( 5.6a ), ( 5.6b| ), and 
( |5.6c ), respectively. Then the divergence of the integrals (5.6) for e > 1 is manifested in the form of the contributions 



divergent as m — > for e > 1. 



As the independent variables in the integrands in (5.6) one can take the mo dulu s k and z = (n • k)/k, the cosine 
between the directions of the momenta p and k. Then the integrals over k in (|5.6| ) in the spherical coordinates are 
written as 



//>oo />00 f-\ 

dK...=S d dKK d - 1 {...)=S d - 1 dK^- 1 dz(l- z 2 ) (d -^l 2 . 

Jo Jo J-l 



(5.8) 



with Sd from Eq. (2.7). Like in Eq. (4.4), the brackets denote the averaging over the unit sphere in d dimensions, 
while the second relation is obtained by the integration over all angles in the (d— l)-dimcnsional subspace orthogonal 
to the vector p. 



Now we turn to the calculation of the coefficient C2 = co\ £ =o in (5.3). The quantity R\ is UV finite and it is sufficient 
to calculate its contribution at e = 0. We put e = in ( 5.6a ) and rewrite the integral in variables z, k. This gives: 



Ri\ e =o = gS d T> = 



gSd- 



16(d-l)(27r)<* J 



dn 



dzil-z 2 )^- 1 ^ 2 



k (d - 1 - 4k 2 + 2ck 2 - 2dzK + 4z 2 n 2 ) 

(1 - 2ZK + K 2 ) d / 2 (l - ZK + K 2 ) 



(5.9) 



The integrals R2 = R3 contain a pole in e (manifestation of the UV divergence in the 1-irreducible function (<^V)i 
see Sec. HI) and can be written as 



R 2 = R 3 = gS d [A/e + 13 + 0(e)]. 



Expression (5.6b) can be written in the form 

gSd-i 



R2 



16(d- l)(27r) rf 7 k 1 



dn 



l+2e 



dz(l- Z 2 ) id - 1)/2 lfj(z,K), 



where 



(1 - d + 3k 2 - dn 2 + 2dzK - 2zk - 2zk 3 ) k 2 
(1 - 2zk + k 2 )(1 - zn + K 2 ) 



(5.10) 



(5.11) 



(5.12) 



The integration region over k can be split in two parts: [0,1] and [l,+oo); the pole in e comes only from the second 
part. The expansion of tp(z, k) at large k has the form ip(z, k) — —2zk + (3 — d — 6z 2 ) + 0(1/ k). The first term is 
odd in z and vanishes after the integration over z; the second term completely determines the residue at the pole: 
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A 



32(d-l)S d 



dz (1 - z 2 )^' 2 (3 - d - 6z 2 ) = ((1 - z 2 )(3 - d - 6z 2 )) 



32(d+2)' 



(5.13) 



where the brackets denote the averaging over the unit sphere in d dimensions and Eqs. (4.4) and (5.6) have been 
used. One can se e th at in representation (5.3), the total pole part 2A/e of the diagrams i? 2 .3 and the pole that comes 
from expression (5.4) can cel e ach other (a consequence of the renormalizability) . Thus the total contribution of the 
diagrams i?i,2,3 into Eq. (5-3) has the form 



c 2 = c 2 | £=0 = 2B + O 



(5.14) 



with B from ( 5.10 ) and T> from (5.9), and for the first two terms of the e expansion of the amplitude R(l, 0, <?*) using 
( |3. 10 ) and (5.3) we obtain 



R{1,0, g*) = (1/2) [1 + 2ac 2 e + 0(e 



(5.15) 



with a = 8(d + 2)/3(d- 1) from (3.10) 



The quantity B can be found as follows. The integrals that remain in (5.11) after the pole part has been subtracted 
(that is, the integral over [0,1] and the integral over [l,+oo) with the substitution ip(z, k) — ► ^{z, k) + 2zk — (3 — d— 6z 2 ) 
in the integrand) are UV finite and we may set e = in them. This gives: 



16(rf-l)S d 
Sd-i 



B 



1 dk 



dk 



— dz(l- z 2 ) {d - 1 ^ 2 i/j(z, k)+ — / dz (1 - z 2 f d ^' 2 yj(z, k) + 2zk - (3 - d - 6z 2 



dk 

T 



dz (1 - z 2 ) {d -^l 2 [ip(z, k) + iP(z,1/k) + 2zn-{3-d- 6z 2 )] . 



(5.16) 



The integrals (5.9) and (5.16) converge and can be evaluated numerically for any given value of d. For d — 3 one 
obtains: 

B = -0.000057, V = 0.06699, c 2 = 0.0669. (5.17) 
The quantities (5.17) for other values of d are given in Table |; we shall discuss them later in Sec. VII. 



VI. FINITE e: NONLOCAL INTERACTIONS, INFRARED SINGULARITIES, AND GALILEAN 

SYMMETRY 



In a number of studies, the results of the RG approach to model (2.1), (2.2), ( |2.9| ) were interpreted and criticized 
in the language traditional for the classical theory of turbulence, comparison with the well-known direct interaction 
approximation (DIA) was made and the validity of the e expansion for finite e ~ 1 was called in question; see e.g. 

§§■ . . „ 

Let us discuss the problem on the example of the pair correlation function in the one-loop approximation (p.q ). 

The corresponding analytical expressions (5.6) involve three vectors p, k, q subject to the restriction p = k + q (here, 
the momenta are not divided by p = |p|). In the following, the external momentum p will be taken to lie in the 
inertial range. Then relevant contributions to ( |5.6| ), roughly speaking, can come from the three regions: k ~ q 3> p, 
p ~ q 3> k, where k is the momentum flowing via the (<pip)o line (any one of the two equivalent ((p<p)o lines for 
the diagram i?i), and p ~ k ~ q (all three momenta lie in the inertial range). We shall refer to the corresponding 
contributions as being determined by UV-nonlocal, IR-nonlocal, and local (in the momentum space) interactions, 
respectively. 

The typical argument against the use of the e expansion for finite e can be formulated as follows. The diagrams 
i?2,3 have poles in e as a manifestation of the UV divergence at e = 0. Thus for small e the leading contribution to 
those diagrams is determined by UV-nonlocal interactions. For e > e c with certain e c = (1) the diagrams become 
IR divergent and do not exist without the IR cutoff m [e c = 1 for all diagrams -Ri.2,3 in ( |5.5| )]. In the "massless" 
model (2.5) with m = they have the form of the poles in e c — e. Physically such IR divergences are explained by 
the so-called sweeping effects, that is, the transport of small eddies as a whole by large ones. This purely kinematic 
phenomenon has no effect on shaping the energy flow over the spectrum and is irrelevant in determining the exponents 
and amplitudes in power laws for Galilean invariant quantities (e.g. equal-time structure functions). According to 
the classical phenomenology the latter are formed by the interactions of the eddies in the inertial range, that is, by 
local interactions (in the above sense); see e.g. B. One may thus believe that the extrapolation of the e expansion 
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(reliable for e <C 1) to finite e ~ 1 overestimates the contribution of the UV-nonlocal interactions and misrepresents, 
or neglects at all, the physically important contributions from the IR-nonlocal and local interactions. 

More speci fical ly, the author of Refs. |ll| extended the well-known direct-interaction approximation pj^|3l|] to 
the problem (2.1), (2.2) with the power-law forcing @.5|). In the DIA, the pair correlation function ((pep) and the 
response function (ipip'} satisfy a closed system of integro-differential equations; their right-hand sides involve the 
skeleton analogs of the diagrams i?i,2,3 with exact ("dressed") lines and bare vertices. The author of Ref. Jlj| used 
a scaling Ansatz, expanded the integrals on the right-hand sides in e and retained their pole parts; this allowed him 
to reproduce some of the results obtained earlier using the RG and the e expansion. On the other hand, one can 
expand the same integrals in the vicinity of the pole in e c — e (for the skeleton one-loop diagrams one has e c = 3/2) 
and obtain a solution whose form is completely different from that of the solution obtained within the e expansion 
p5| . Its extrapolation to the physical value reproduces the well-known solution [^9| of the original DIA equations. 
Although the e expansion fails to reproduce this solution, no definitive conclusion can be made at this point: for the 
different-time pair correlation function, neither of the aforementioned solutions is correct for e > s c . 

Indeed, it has long been realized that the DIA misrepresents the contributions of the IR-nonlocal interactions, at 
least in the energy spectrum (it is determined by the equal-time pair correlation function and is Galilean invariant); see 
Ref. [p9|-[3l|| . Although the original Dyson- Wyld equations for the correlation and response functions are exact, they 
involve infinite series of skeleton diagrams; the practical solution implies a truncation (the one-loop truncation gives 
the DIA) followed by the substitution of a scaling Ansatz. As a result of the approximation, the Galilean symmetry 
is violated, and the IR singularities, related to the sweeping effects, do not cancel out in Galilean invariant quantities 
like the energy spectrum]]] This leads to spurious strong dependence of the latter on the IR scale and erroneous value 
(3/2 instead of 5/3) for the corresponding exponent. No closed equation, however, can be written for the equal-time 
correlation function, and the problem cannot be solved by taking into account more diagrams. 

The key difference between the RG and the self-consistency approaches is that the former is based on the ordinary 
perturbation t heor y, tha t is , ex pansion in the nonlinearity in fl2.1| ), which involves a formal expansion parameter go. 
The problem (2.1), ( |2.2| ), ( [2~9| ) is Galilean invariant for any go, e and d, so that the perturbation theory in go is 
manifestly Galilean invariant in any given order for any value of e and d; this is equally true for the renormalized 
perturbation theory (in g). The use of the RG and OPE implies some infinite resummations of the original pertur- 
bation series within controlled approximations and therefore does not violate the Galilean symmetry of the original 
perturbation expansion. 

The plain e expansion is indeed not suitable for the description of the IR singularities related to the sweeping. 
Even in the standard cj> 4 model of critical behavior, where the IR singularities are weaker, e expansions of the form 
( |3.16 ) cannot be used for the analysis of the asymptotical behavior at m — > 0. As already discussed in the end of 
Sec. Ill , this fact does not hinder the use of the RG techniqu e, wh ich should be combined with the operator product 
expansion to derive resummed representations of the form ( [3.171 ); see e.g. The distinguishing feature of our 

model (2.1), (2.2), (2.9) is the existence in Eq. (3.17) of negative dimensions Ap < 0. The summation of the most 
singular contributions coming from the operators ip n was performed in Ref. p5| , the generalization to the case of a 
time-dependent large-scale field is given in p8|] . This gives the adequate description of the sweeping effects within the 
RG formalism; see also Refs. P,fTo|]. 

Admittedly, negative dimensions can be reliably identified, and their contributions can be summed up, only with the 
aid of nonperturbative methods that allow one to go beyond the plain e expansion: functional Schwinger equations, 
Ward identities that express Galilean symmetry, and infrared perturbation theory. These techniques are discussed in 
Refs. |],[l(| in detail. At the moment, however, we are interested in the equal-time pair correlation function (|5.l|). 
As already mentioned in the end of Sec. [II, the operators ip n , as well as other noninvariant operators, give no 
contribution to the corresponding representation (3.17) in agreement with the fact that Galilean invar iant quantities 
are not affected by the sweeping. Only the dimensions of Galilean invariant operators can appear in ( |3.17 ), and the 
singular behavior can be related to invariant operators with < 0. No such operator, however, has been presented 
in any study we know of, and for the time bei ng w e can assume that all dimensions Ap are non-negative. This fact 
suggests that the scaling function i?(l,it,g*) (5^2) in Sec. |V| is finite at u ~ and the value of i?(l,0,<?*) can be 
calculated within the framework of the "massless" model ( |2.5| ) and the e expans ion. 

The cancellation of the IR divergences for 1 < e < 2 in the sum of diagrams ( |5.6|) was demonstrated in Sec. fv| fo r 
general d. Below we consider them in the limit d — * oo, where explicit expressions for the integrals entering Eq. (5.6) 



1 The DIA is not a given-order approximation to the exact Dyson- Wyld equations in any small, at least formal, expansion 
parameter. It can be understood as the leading approximation for N — > oo in certain extension of the original NS equation to 
the case of N interacting velocity fields J32|, but such an extension is not Galilean invariant unless N — 1. 
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can be written for all < s < 2 (in Sec. [y], only two terms of the e expansion were calculated). 

In the diagrams i?2 = R3 nontrivial dependence on d comes only from the angular averaging. It follows from Eqs. 
(4.4) that for d — > 00, the angular averages behave as (z 2n ) oc d~ n . Thus in order to find the leading term of the 
large-d behavi or, it is sufficient to discard all such averages with n > 1 , or, equivalently, to put z = (n • k)/fc = in 
the integrand (5.6b). In the variables k, z introduced above Eq. (5.8) this leads to the replacements (1 — z 2 ) — > 1, 



(1 + ft 2 ), 1 + k 2 + q 2 — > 2(1 + k 2 ), which along with the asymptotic relation Q 2 — — d(l + n 2 ) + 0(d°) for the 

-gS d f°° J ^ 2e -g§. 



quantity Q2 in (5.7) gives 
R2 + R3 — 2i?2 



-9 
8(2tt) 



„2-d-2e 



dn 



dn 



(1 



16 s 



T(l + e)r(l-e). 



(6.1) 



The latter equality in ( S.l) holds for e < 1; for e > 1 the integral over k diverges at small k. 

In the integrand ( 5.6a ) for diagram nontrivial dependence on d come s also from the factor q~ d and the procedure 
described above cannot be applied. In the variables k, z expression (5.6a) takes on the form 



Ri 



gSd- 



8(d-l)(2nY J 



dK / dz(l-z 2 )^' 2 



Q 1 K 1 - 2e q- d - 2e 
(1 + k 2 + q 2 ) 



(6.2) 



with q 2 = 1 



(1 - z 2 ) d ' 2 q- d 



2zk and Q\ from (5.7). For large d one can write 

, {k-z) 2 



(1- Z 2 ) d / 2 (1+K 2 



2zk) 



-d/2 



(1-z 2 ) 



-d/2 



exp 



d (k — z) 2 
2 (1-z 2 ) 



Along with the relat ion Sd-i/Sd — d/2it [1 + 0(1/ d)], which follows from the Stirling formula for the T function 
that enters Eq. (2.7), this gives: 



(1 _ jy.d-D/iq-dg^ „ S d 5(K - z). 



We substitute this relat ion into (3.2), retain only intervals [0, 1] in the integrals, use the relation Qi\k= z — d + 0(dP) 
for the quantity Q\ in (5.7) and perform the integration over any one of the variables; this gives: 



Ri 



gSd 



8(d-i) 



1 ,1 

dn / dzQi 
a Jo 



K l-2e q -2e gSd 
(l + K 2 + q 2 ) [K Z) ~ 16 



dn 



k 1 - 2 ' = -gS d T 2 (l-e) 
(1-k 2 Y 32 T(2-2e)' 



(6.3) 



Like for the integral (3.1), the latter equality holds for e < 1; for e > 1 the integral over k diverges at small k — > 
and K — * 1. Adding together the expressions (|SJ|), ( |6.3| ) with the loopless contribution (f>A\), where the limit d — » 00 
of the explicit form (3.6) should be used for the coefficient an, gives: 



R(l,0,g) 



gSp 

32 



11(e), 11(e) 



2 2, N , . T 2 (l-e) 



(6.4) 



(the product gSd at the fixed point has a finite limit for d = 00; see Eq. ( 3. 10f )) . The sum of the integrals ( p.lD , ( |6.3| ) 
converges for all < e < 2: the divergences in t he integrals ( |6.l[ ) and (6.3|) at K — > cancel each other, while the 
cancellation of the divergence in the integral ( |6.3|) a t k — > 1 with the divergence in the diagram i? 3 becomes obvious 
after the replacement /c — > 1 — k in expression ( p.l| ) for i?2- Thus the expression ( |6.4| ) for the sum of the integrals 
holds for the whole interval < s < 2. 

In Eq. JS-4); the cancellation of the poles at e = (a consequence of the renormalizability) and e = 1 (a consequence 
of the Galilean symmetry) are obvious; the renormalizability and the Galilean invariance of the model guarantee similar 
cancellations in all the higher orders of the (renormalized) perturbation theory. 

In models of quantum field theory, the elimination of UV divergences is needed to ensure existence of meaningful 
(finite) results at the physical value e — 0. In our case it is not necessary to take the limit e — * and ensure that it 
exists, that is, there is no need to renormalize the model. However, it turns out that the possibility itself of doing 
this implicitly contains valuable i nfor mation about the original model with finit e e: it is this possibility that ensure 
the validity of the RG equations (3/7) and the operator product expansion ( 3.17 ). 

Physically, elimination of the UV singularities (poles in e) removes the contributions of the UV-nonlocal interactions 
that would obscure relevant local interactions. Coefficients of the renormalized perturbation series for the correlation 
functions have a finite limit at e — *■ 0: they are determined by the local inte ract ions, while the UV-nonlocal interactions 
are taken into account in the formulas of multiplicative renormalization (3.4). A fine correlation that exists between 
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different orders of the perturbation theory is expressed by the RG equations (3.7); solving them allows one to perform 
certain infinite res ummation of the perturbation series for (renormalized) correlation functions. 

Expression ( |6.4[ ) contains a pole at e — > 2: R(s) — 6/(2 — e) [for general d we would obtain R(e) ~ 2(3<i 2 + d — 
12) /d(d + 2) (2 — e)]; higher-order poles will necessarily appear in the higher orders of the perturbation theory. Such 
singularities are necessary to ensure the exi stence of the finite limit for the pair correlation function at s — > 2 in 
terms of the physical parameter £ from ( |2.4[) . For the correlation function of n fields ip, they must combine into the 

s ing ularity (2 — e)~"/ 3 that will cancel the vanishing factor D^ 3 oc £™ 3 (2 — e)"/ 3 in representation ( |3 .15 ); see Eq 



fl2.8| ). The conjecture that the singularities at e — > 2 indeed combine into the needed expressions is confirmed by 
two examples: the pair correlation function in the exactly solvable Heisenberg model |E6j and the triple correlation 



function in our case; see the exact expression (7.11) in Sec. VII 



VII. TWO-LOOP CALCULATION OF THE KOLMOGOROV CONSTANT 



The Kolmogorov constant Ck can be defined as the (dimensionless) coefficient in the inertial-range expression 
S^M = Ck (£r) 2 / 3 for the second-order structure function, predicted by the Kolmogorov-Obukhov theory and con- 
firmed by experiment |l]||- Here £ is the mean energy dissipation rate and the n-th order (longitudinal, equal-time) 
structure function is defined as 

S n (r) = ([<p r (t,x + r)-(p. r (t,x)] n ), <Pr = (v>i'n)/r, r = \r\. (7.1) 
Alternatively, the Kolmogorov constant C' K can be introduced via t he re lation E(k) — C' K £ fc -5 / 3 , where the 



energy spectrum E(k) is related to the equal-time correlation function (5.1) as E(k) — Sd(d — l)k d 1 G(k)/2. From 



the definitions one can derive the following relation between these two constants: 



which for d — 3 gives (cf. W 



_ 3-2V 3 r(2/3)r(d/2) 
CK " (d + 2/3)r(d/2 + l/3) 0/0 [J - ' 



_ 9-2V 3 y^r(2/3) _27 
Ck ~ 22T(ll/6) ° K ~ 55 T{1/3)Ck - (? - 3) 



Using the exact relation S3 (r) = — I2£r/d(d + 2) that follows from the energy balance equation (see e.g. for 
d = 3), the constant Cjf can be related to the inertial-range skewness factor S: 

S = Sz/Sl 12 = -[l2/d{d + 2)]C K 3/2 . (7.4) 



All these relations refer to real physical quantities in the inertial range, which in the stochastic model (|2.l|), (2.2) 
corresponds to e = 2 and m = in the random force correlator 

Many studies have been devoted to the derivation of Ck within the framework of the RG approach; see Refs. 
P, ^0|j2qj3^ -p9t. In order to obtain Ck, it is necessary to augment the solution of the RG equation for S2 by some 



formula that relates the amplitude Do in the random force correlator (2.9) to the physical parameter £. In particular, 
in |3^j34| the first-order term of the e expansion for the pair correlator was combined with the so-called eddy-damped 
quasinormal Markovian approximation for the energy transfer function, taken at e = 2. More elementary derivation, 



based on the exact relation (|2.4[ ) between £ and the function df{k) from (2.9) was given in p5[ ; see also §,0. In 



spite of the reasonable agreement with the experiment, such derivations are not immaculate from the theoretical 
view poin ts. Their common flaw is that the relation between £ and Dq is unambiguous only in the limit e — * 2 [see 



Eq. fl2.8|) in Sec. [Ill, s o that the coefficients of the corresponding e expansions can, in fact, be made arbitrary; see the 
discussion in Ref. J37| and Sec. 2.10 of |l(J. The ambiguity is a consequence of the fact that the notion itself of the 
Kolmogorov constant has no unique extension to the nonphysical range e < 2. 

The experience on the RG theory of critical behavior shows that well-defined e expansions can be written for 
universal quantities, such as critical exponents, normalized scaling functions and ratios of amplitudes in scaling laws 
The constant Ck extended to the range e < 2 as in |53|,[54| or [p5| involves a bare parameter, Dq, and hence is 
not universal in the above sense. 

To circumvent this difficulty, we propose below an alternative derivation that relates Ck to an universal quantity 
and thus does not involve any relation between Dq and £. Consider the ratio 
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Q{e) = V r S 2 {r)/\S 3 {r)\ 2 ' 3 = V r S 2 (r)/(-S 3 (r)) 2 / 3 



(7.5) 



with T> r = rd/dr. As we shall see below, the quantity Q(e) is universal and can be calculated in the form of a 
well-defined e expansion. On the other hand, its value at e — 2 determines the Kolmogorov constant and the skewncss 
factor through the exact relations 



C K = [3Q(2)/2] [12/d(d+2)] 



2/3 



S = - [3Q(2)/2]' 



-3/2 



that follow from the definitions, relation ( |7.4| ) and the identity V r r a = ar a for any a. 

Solving the RG equat ion s for the quantities in Q(e) in the inertial range (Ar — > oo and mr 
< e < 2 gives (see Sec. pl|): 



(7.6) 



0) for general 



S 3 (r) = D r- 3A -f 3 (e), V r S 2 {r) = D 



2/3„-2A 




f/a(e), A„ = l-2e/3. 



(7.7) 



for 



It is important here that the operation T> r kills the constant contribution (ip 2 ) in 62 that diverges as A 
e < 3/2 (see below); i n S3 such constant contributions are absent. 

It follows from Eq. (7.7) that the ratio Q(e) = fz/{— f 3 ) 2 ^ 3 , in contrast to its numerator and denominator, does not 
depend on the amplitude Dq, it is universal (in the above sense) and can be calculated in the form of a well-defined 
e expansion. The latter has the form Q(e) = e 1 ^ 3 p(e), where p(e) is a power series in s (s ee be low). We s hall find 



the first two terms of p{e), which corresponds to the two- loop accuracy in representations ( 4.23 ) and ( 5.15 ). In this 
sense, one can speak about the two-loop approximation for the Kolmogorov constant (previous attempts have been 
confined with the first order). 

To avoid possible misunderstandings, we stress again that we shall not try to extend the definition of the physical 
quantities Ck and S to the whole interval < e < 2 and to construct the corresponding e expansions from the known 
expansion for Q(e). I nste ad, the latter is used to give the value of Q(2), which, in its turn, determines Ck and S 
through the relations (7.6) that make sense only for e = 2. 



Using the definition (7.1), the function S 2 can be related to the momentum-space pair correlation function (5.1) as 
follows: 



S 2 (r) = 2 J G{k) [1 - (k • r) 2 /(fcr) 2 ] {1 - exp [i(k • r)]} . 

Applying the operation T> r = rd/dr gives: 

V r S 2 (r) = 2 J G(k) [1 - (k • r) 2 /(fcr) 2 ] (k • r) sin(k • r). 



(7.8) 



(7.9) 



In order to obtain the inertial-range form of the function T> r S 2 (r) it is sufficient to substitute the asymptotic expression 
(5.2) into (7.9). A straightforward calculation gives: 



V r S 2 (r) 



2(d - l)r(2 - 2g/3) 
(47r) d / 2 r(d/2 + 2e/3) 



gl /3 R(l,0,9*)D 2 /3 (r/2Y 



2A« 



(7.10) 



with the amplitude i?(l,0,g*) from (5^2). It is important here tha t th e resulting integral exist s for all < e < 2. It 
i s th e operation T> r that ensures the convergence of the integral (7.9) with the function (5.2); the original integral 
(7.8) would be UV divergent for < e < 3/2 while the analogous expression for the pair correlator (that is, the plain 
Fourier transform of ([5.2|)) would be IR divergent for e > 3/2. 

The needed terms of the e expansion for f 3 can be obtained not only from the direct perturbative calculation, but 
also from the exact expression 



3(d - l)r(2 - e) 
(47r) d / 2 r(d/2 + e) 



A>(r/2) 



-3A„ 



(7.11) 



that follows from the energy balance equation and in the limit e — > 2, along with the formula ( |2.8| ), reproduces the 
correct coefficient —12/d(d+ 2). 

The pole ~ 1/(2 — e) cancels the vanishing factor ~ (2 — s) in_(2.8) and ensures the finiteness of S3 at e — 2 



when expressed in terms of the physical parameter £. In the ratio (7.5) this pole must cancel with the (hypothesized) 
singularity in the p air co rrelat or; se e the discussion in the end of Sec. VI. 
Thus from Eqs. (7.10) and (7.11) for the ratio Q(e) one obtains 
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Q(e) = [4(d - l)gJ d /9] 1/3 A(e, d)R(l, 0, <?*), 



(7.12) 



where the coefficient 



Me d) = r(2-2 e /3)rV3( d/ 2)r'/3( d/2 + e ) _ 

A(£ ' d > - T(d/2 + 2e/3)r2/3(2 -e) + ° {£ ' 



(7.13) 



has no term of order 0(e) and in our approximation can be replaced by unity. Substituting the two-loop result ( 4.23| ) 
for t he co ordinate of the fixed point g„ and the first two terms (5.15) of the e expansion of the amplitude R(l, 0, g*) 
into ( 7.12 ) gives 



Q(e) = (l/3)[4e(d+2)] 1 / 3 [l + e (A/3 + 2ac 2 ) + 0(e 2 )} 



(7.14) 



with a. and A from Eq. ( [1.23] ) and c 2 from ( 5.14 ). For d = 3 this gives 

Q(e) = (l/3)(20e) 1/3 [l + 0.525e + 0(e 2 )] . (7.15) 
The value of Q(e) at e = 2 determines the Kolmogorov constant and skewness factor through the exact relations 



(Q, which for d = 3 have the form C K = 6 • 10~ 2 / 3 Q(2) and 5 = -[1.5 ■ Q(2)]~ 3 / 2 . Let us denote as Cg' and 
the results obtained using the n-loop approxima tion for the RG functions and the corresponding approximations for 
the scaling functions. Substituting the value of ( 7.15 ) with e = 2 into these relations we obtain 



C 



3.02, 5 (2) 



-0.15. 



If we had retained only the first-order term in (7.15) we would have obtained 



C%> = 1.47, S {1) = -0.45. 
We also recall the experimental estimates recommended in [jjj : 

C K ~ 1-9, S w -0.28; 

it is worth noting that they lie in between of the first-order and second-order approximations. 

Let us conclude this section with a brief discussion of the d-dimensional case. Results of the two-loop calculation 
for several values of <i, including the limits d — > oo and d — > 2, are given in Table |j| It includes: the parameter A that 
determines the second-order corrections to the coordinate of the fixed point g* from Eq. ( 4.23| ) and the exponent ui 
from Eq. (4.25), the parameter C2=_^\ e =o that determines the second-order correction to the scaling function of the 
equal-time pair correlator in Eq. ( [5.15| ), the contributions B and T> from the diagrams i?2,3 and R\ that determine c 2 
via Eq. ( 5.14 ), and the first-order and second-order approximations for the Kolmogorov constant (n = 1,2). 

One can see that, for d — > oo, the two-loop contributions become relatively smaller for all these quantities. This 
gives some quantitative support to the old idea that 1/d can be a promising expansion parameter for the stochastic 
NS problem pQ ]. However, in contrast to the Kraichnan model, where the 0(l/d) approximation for the anomalous 
exponents is a function linear in e fl9|| , in our case the second-order relative corrections have a finite limit at d — > oo, so 
that the 0(1/ d) approximation for g*, u, Q(e) are infinite series in e. On the other hand, like for the Kraichnan model, 
the calculation of the diagrams drastically simplifies for d — > oo: the leading terms can be calculated analytically. We 
believe that the three-loop calculation is a feasible task at d — oo; this work is now in progress. 

For d — > oo, the constant Ck decreases as Ck oc 1/d; see Table § From Eq. one then obtains C' K oc d 1 ^ 3 for 

the coefficient in the energy spectrum, in agreement with the earlier results obtained within the DIA |40| or the RG 

For d — > 2, the two-loop contributions diverge; this is a manifestation of the additional UV divergence that emerge in 
the 1-irreduciblc function (tp'ip 1 ) at d = 2. In this region our results cannot be trusted, and additional renormalization 
should be performed in order to remove the divergences at d = 2; see e.g. Sec. 3.10 of Ref. [[[o| and references therein. 
This is necessary, in particular, in the discussion of possible crossovers in the inertial-range behavior that can occur 
bet ween d = 3 and 2; see Refs. 41 . It is worth noting that, for all d, the main contribution to the coefficient c 2 = V+2B 
in ( |5.14 ) comes from the term V, that is, from the diagram R\. What is more, the contribution from B totally vanishes 
at d = oo while the contribution from T> diverges at d = 2. It thus might happen that the aforementioned additional 
renormalization is needed to obtain more accurate numerical predictions in three dimensions; we shall consider this 
important problem elsewhere. 
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VIII. CONCLUSION 



We have accomplished th e com plete two-loop calculation of the renormalization constant and RG functions for 
the stochastic problem (2.1 )— ( 2.9) and derived the coordinate of the fixed point, the UV correction exponent u>, 
the Kolmogorov constant Ck and the inertial-range skewness factor S to the second order of the corresponding e 
expansions. The new point is not only the inclusion of the second-order corre ctio n, but also the derivation of Ck 
through the universal (in the sense of the theory of critical behavior) quantity (7.5). 

Of course, one should have not expected that the second-order terms of the s expansions would be small in 
comparison to the first-order terms. The experience on the RG theory of critical behavior shows that such corrections 
are not small for dynamical models (in contrast to static ones) and for amplitudes (in contrast to exponents); see ||,[l). 
It is thus rather surprising that in our case the account of the two-loop contributions leads to reasonable changes in 
the results. 

Although the e 2 correction to uj in ( 1.25 ) is rather large, it does not change its sign and hence does not destroy the 
IR stability of the fixed point for the real e — 2. 

The first-order approximation Cjp = 1.47 underestimates, and the second-order approximation — 3.02 overes- 
timates the conventional experimental value of the Kolmogorov constant Ck ~ 1-9 Jlj]. Thus the experimental value 
of Ck (and hence for S) lies in between of the two consecutive approximations. A similar situation is encountered for 
the well-known Heisenberg model Q , where the analog of the Kolmogorov constant is known exactly and lies between 
the first-order and second-order approximations given by the corresponding e expansion If we assume, by the 
analogy with the Heisenberg model, that the (unknown) exact predictions for Ck and S lie between the first two 
approximations, we may conclude that our calculation has given satisfactory estimates for these quantities. 

The two-loop corrections become relatively small for d — > oo, which confirms the relevance of the 1/d expansion for 
the issue of inertial-range scaling behavior. 

One may thus conclude that our results confirm the applicability of the RG formalism and the e expansion to the 
calculation of important characteristics of the real fluid turbulence. The next important steps should be the derivation 
of the two-loop results for the case of a passive scalar advected by stochastic NS equation; inclusion of the additional 
renormalization near two dimensions, and the three-loop calculation in the limit d — > oo; this work is now in progress. 
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APPENDIX: 



Below we give explicit expressions for the integrands /0-/4 in the expressions ( |4.8| ) for the two-loop diagrams of the 
1-irreducible correlation function (ip'<p). They have the form 

/,(;,,) = ,= (1 - u{d _^ d + 2)Sd [AM + *(-.,.)] 



with Sd from ( |2.7j ) and x = k\/k2- Note that the functions fi(z,x) are even in z [see the remark below Eq. (4.8)] 
The functions tpi = ipi(z, x) in d dimensions have the forms (the function ipg is even in z, while the other functions tp. 
are not) 

^ = {(x 20 + l)[(-3d 3 + 8d 2 - Id- 16) + (-14d 2 + 100d + 40)z 2 - \2Mz 4 ] 
+ (x 18 + x 2 ) [(-30d 3 + 76d 2 - 78rf - 64) + (12d 3 + 24rf 2 + 552d - 144)z 2 
+ (-184d 2 - 76<i+ 184)z 4 - UOdz 6 ] 

+ (x 16 + x 4 ) [(-135d 3 + 328d 2 - 395d + 48) + (96d 3 + 522rf 2 + 1080(2 - 1272) 
+ (63d 3 - 1796d 2 + 1495d + 1728)z 4 + (1054<i 2 - 2984<i - 576)z 6 + 816dz 8 ] 
+ (x 14 + x e ) [(-360d 3 + 848d 2 - 1160d + 768) + (336rf 3 + 1952d 2 - 88d - 2752) 



z 2 



z 2 
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+ (378d 3 - 7348d 2 + 5918d + 3136)z 4 + (-312rf 3 + 6124rf 2 - 5320rf - 1232) 2 6 
+ (-1432d 2 - 208d+ 224)z 8 + 1248rfz 10 ] 

+ (x 12 + x e ) [(-630d 3 + 1456rf 2 - 2158rf + 2016) + (672d 3 + 3844rf 2 - 3740rf - 2864)z 2 

+ (945d 3 - 16172d 2 + 10945d - 64)z 4 + (-1248rf 3 + 15362rf 2 - 2296d + 1408)z 6 

+ (240d 3 - 5376d 2 - 4224d- 768)z 8 + (832d 2 + 1792rf + 128)z 10 - 640dz 12 ] 

+ 2a; 10 [(-378d 3 + 868d 2 - 1322rf + 1344) + (420rf 3 + 2376d 2 - 3024d - 1200)z 2 

+ (630rf 3 - 10436d 2 + 6478rf - 1656)z 4 + (-936d 3 + 10292rf 2 + 416d + 2064)z 6 

+ (240rf 3 - 3912d 2 - 3360d - 480)z 8 + (832d 2 + 544rf)z 10 + (-128d 2 + 256d)z 12 ] } 

/ [(x 2 + 2xz + l) 3 (x 2 + xz + l) 3 (x 2 - 2xz + l) 3 (x 2 -xz + l) 3 ] ; 

ip! = 2{x 7 z(-2d 2 + 2d) + x 6 [(d 3 -Ad 2 + d + 6) + (-6d 2 + 6d - I2)z 2 } + x 5 z[(5d 3 - 20d 2 + 7d - 4) 
+ (-Ad 2 + Ad- 16)z 2 ] + x 4 [{3d 3 - I0d 2 + 3d) + (8d 3 - 26d 2 + lOd - 28)z 2 ] 
+ x 3 z [{10d 3 - 28d 2 + lOd - 20) + (4rf 3 - 8rf 2 + Ad)z 2 } + x 2 [{3d 3 - 8d 2 + 3d - 6) 
+ (8d 3 - 16rf 2 + 8d)z 2 ] + xz(5d 3 - 10d 2 + 5d) + d 3 - 2d 2 + d) 
I [(z 2 + 2a ; z + l) 2 (:E 2 + .TZ + l) 2 ] ; 

ip 2 = 2{x i (2d 2 -2d + 2)+ x 3 z(-4d 2 + U - 4) + x 2 [(5d 2 - 5d + 6) + (2d 2 -2d+ 2)z 2 ] 
+ xz(-hd 2 + bd - 6) + 3d 2 - 3d + 6} [2x 3 z + x 2 (d-3) + xz(-2d + 2) + d - l] 
/ [(x 2 - 2xz + l)(x 2 -xz + l) 3 ] ; 

iP 3 = 2d(d - 1) [2x 3 z + x 2 (d-3) + xz(-2d +2)+d- l]/[(x 2 - 2xz + l)(x 2 - xz + 1)] ; 

V> 4 = { (x 4 + 1) (3d 2 -3d +2) + (x 3 + x)z(-7d 2 + 7d-8) + 2x 2 [(3d 2 -3d +2) + (2d 2 -2d + 4)z 2 ] } 
[ (x 4 + l)(-d + 1) + (x 3 + x)z(2d - 4) + 2x 2 (-d + 3)]/ [(x 2 - 2xz + l) 3 (a; 2 - xz + l) 3 ] . 



[1] A. S. Monin and A. M. Yaglom, Statistical Fluid Mechanics (MIT Press, Cambridge, Mass., 1975), Vol.2. 
[2] U. Frisch, Turbulence: The Legacy of A. N. Kolmogorov (Cambridge University Press, Cambridge, 1995). 
[3] J. Zinn- Justin, Quantum Field Theory and Critical Phenomena (Clarendon, Oxford, 1989). 

[4] A. N. Vasil'ev, Quantum-Field Renormalization Group in the Theory of Critical Phenomena and Stochastic Dynamics 
(St Petersburg Institute of Nuclear Physics, St Petersburg, 1998) [in Russian; English translation: Gordon and Breach, in 
preparation]. 

[5] D. Forster, D. R. Nelson, and M. J. Stephen, Phys. Rev. Lett. 36, 867 (1976); Phys. Rev. A 16, 732 (1977). 
[6] C. De Dominicis and P. C. Martin, Phys. Rev. A 19, 419 (1979); Progr. Theor. Phys. Suppl. No 64, 108 (1978). 
[7] J. D. Fournier and U. Frisch, Phys. Rev. A 19, 1000 (1983). 

[8] L. Ts. Adzhemyan, A. N. Vasil'ev, and Yu. M. Pis'mak, Theor. Math. Phys. 57, 1131 (1983). 

[9] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil'ev, Usp. Fiz. Nauk, 166, 1257 (1996) [Phys. Usp. 39, 1193 (1996)]. 
[10] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasiliev, The Field Theoretic Renormalization Group in Fully Developed 

Turbulence (Gordon & Breach, London, 1999). 
[11] R. H. Kraichnan, Phys. Fluids A 30, 2400 (1987). 
[12] S. Chen and R. H. Kraichnan, Phys. Fluids A 1, 2019 (1989). 

[13] E. V. Teodorovich, Izv. Akad. Nauk SSSR, Ser. Mek. Zhidk. Gaza, No 4, 29 (1987) [in Russian]; Sov. Phys. Doklady 33, 
247 (1988); Izv. Akad. Nauk, Fiz. Atmosf. Okeana, 29, 149 (1993) [in Russian]. 

[14] E. V. Teodorovich, On the infrared divergences and the role of local and nonlocal interactions in forming developed turbu- 
lence. Preprint No. 388, IPM, USSR Academy of Science (Moscow, 1989) [in Russian]. 

[15] S. L. Woodruff, Phys. Fluids A 4(5), 1077 (1992); Phys. Fluids A 6(9), 3051 (1994). 

[16] R. H. Kraichnan, Phys. Fluids 11, 945 (1968). 

[17] R. H. Kraichnan, Phys. Rev. Lett. 72, 1016 (1994); 78, 4922 (1997). 

[18] K. Gawedzki and A. Kupiainen, Phys. Rev. Lett. 75, 3834 (1995); D. Bernard, K. Gawedzki, and A. Kupiainen, Phys. 
Rev. E 54, 2564 (1996). 



19 



[19] M. Chertkov, G. Falkovich, I. Kolokolov, and V. Lebedev, Phys. Rev. E 52, 4924 (1995); M. Chertkov and G. Falkovich, 

Phys. Rev. Lett. 76, 2706 (1996). 
[20] U. Frisch, A. Mazzino, and M. Vergassola, Phys. Rev. Lett. 80, 5532 (1998); Phys. Chem. Earth B 24, 945 (1999); U. Frisch, 

A. Mazzino, A. Noullez, and M. Vergassola, Phys. Fluids 11, 2178 (1999). 
[21] G. Falkovich, K. Gawedzki, and M. Vergassola, Rev. Mod. Phys. 73, 913 (2001). 
[22] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil'ev, Phys. Rev. E 58, 1823 (1998). 

[23] L. Ts. Adzhemyan, N. V. Antonov, V. A. Barinov, Yu. S. Kabrits, and A. N. Vasil'ev, Phys. Rev. E 63, 025303(R) (2001); 

64, 019901(E) (2001); E 64, 056306 (2001). 
[24] N. V. Antonov, A. Lanotte, and A. Mazzino, Phys. Rev. E 61, 6586 (2000); L. Ts. Adzhemyan, N. V. Antonov, and 

A. V. Runov, Phys. Rev. E 64, 046310 (2001) 
[25] L. Ts. Adzhemyan, N. V. Antonov, and A. N. Vasil'ev, Zh. Eksp. Teor. Fiz. 95, 1272 (1989) [Sov. Phys. JETP 68, 733 

(1989)]; N. V. Antonov, Zapiski Nauchnykh Seminarov LOMI 169, 18 (1988) [English translation: Sov. Math. Zapiski 

LOMI]. 

[26] L.Ts. Adzhemyan and N.V. Antonov, Theor. Math. Phys. 115, 562 (1998). 

[27] L. Ts. Adzhemyan, A. N. Vasil'ev, and M. Hnatich, Theor. Math. Phys. 74, 115 (1988); N. V. Antonov, S. V. Borisenok, 
and V. I. Girina, Theor. Math. Phys. 106, 75 (1996); N. V. Antonov and A. N. Vasil'ev, Theor. Math. Phys. 110, 97 
(1997). 

[28] L. Ts. Adzhemyan, N. V. Antonov, and T. L. Kim, Theor. Math. Phys. 100, 1086 (1994). 
[29] R. H. Kraichnan, J. Fluid. Mcch. 5, 497 (1959). 

[30] S. A. Orszag, in Fluid Dynamics, Les Houches, 1913 (Eds. R. Balian, J.-L. Peube) (London: Gordon and Breach, 1977) 
p.235. 

[31] W. D. McComb The Physics of Fluid Turbulence (Oxford: Clarendon, 1990). 

[32] C.-Yu. Mou and P. B. Weichman, Phys. Rev. Lett. 70, 1101 (1993); Phys. Rev. E 52, 3738 (1995). 

[33] V. Yakhot and S. A. Orszag, Phys. Rev. Lett. 57, 1722 (1986); Journ. Sci. Comp. 1, 3 (1986). 

[34] W. P. Dannevik, V. Yakhot, and S. A. Orszag, Phys. Fluids 30, 2021 (1987). 

[35] Y. Zhou, G. Vahala, and M. Hossain, Phys. Rev. A 40, 5865 (1989). 

[36] D. Carati, Phys. Rev. A 41, 3129 (1990); Phys. Fluids A 2, 1854 (1990). 

[37] S. H. Lam, Phys. Fluids A 4(5), 1007 (1992). 

[38] M. J. Giles, Phys. Fluids 6(2), 595 (1994). 

[39] Y. Nagano and Y. Itazu, Fluid Dynamics Research, 20, 157 (1997). 

[40] U. Frisch, J. D. Fournier, and H. A. Rose, J. Phys. A: Math. Gen. 11, 187 (1978). 

[41] U. Frisch, M. Lesieur, and P. L. Sulem, Phys. Rev. Lett. 37, 895 (1976); J. D. Fournier and U. Frisch, Phys. Rev. A 17, 
747 (1978); V. Yakhot, Phys. Rev. E 63, 026307 (2001). 



TABLE I. Results of the two-loop calculation for different values of d. 



d= 2 + 2(5 d = 2.5 d = 3 d = 5 d -> oo 

~ -1/35 + 0(5") -2.296 -1.101 -0.560 -l/3 + 0(l/d) 

B O(5 ) 0.0013 -0.000057 -0.00194 0(l/d) 

V 1/645 + 0(6 ) 0.0999 0.06699 0.0436 l/32 + 0(l/d) 

c 2 1/645 + 0(5°) 0.103 0.0669 0.0397 l/32 + 0(l/d) 

C^P 2.08 + 0(5) 1.72 1.47 1.35 5.24/d + 0(l/d 2 ) 

0.93/5 + 0(5°) 4.74 3.02 1.84 5.82/d + 0(l/d 2 ) 
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